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Abstract 

Molecular phylogenetic techniques do not generally account for such common evolutionary events as site 
insertions and deletions (known as indels). Instead tree building algorithms and ancestral state inference proce- 
dures typically rely on substitution-only models of sequence evolution. In practice these methods are extended 
beyond this simplified setting with the use of heuristics that produce global alignments of the input sequences — 
an important problem which has no rigorous model-based solution. In this paper we open a new direction on 
this topic by considering a version of the multiple sequence alignment in the context of stochastic indel models. 
More precisely, we introduce the following trace reconstruction problem on a tree (TRPT): a binary sequence 
is broadcast through a tree channel where we allow substitutions, deletions, and insertions; we seek to recon- 
struct the original sequence from the sequences received at the leaves of the tree. We give a recursive procedure 
for this problem with strong reconstruction guarantees at low mutation rates, providing also an alignment of 
the sequences at the leaves of the tree. The TRPT problem without indels has been studied in previous work 
(Mossel 2004, Daskalakis et al. 2006) as a bootstrapping step towards obtaining information-theoretically opti- 
mal phylogenetic reconstruction methods. The present work sets up a framework for extending these works to 
evolutionary models with indels. 

In the TRPT problem we begin with a random sequence xi,...,Xk at the root of a d-ary tree. If vertex 
v has the sequence j/i, . . . yu v , then each one of its d children will have a sequence which is generated from 
yi , . . . yk v by flipping three biased coins for each bit. The first coin has probability p s for Heads, and determines 
whether this bit will be substituted or not. The second coin has probability pd, and determines whether this bit 
will be deleted, and the third coin has probability pi and determines whether a new random bit will be inserted. 
The input to the procedure is the sequences of the n leaves of the tree, as well as the tree structure (but not 
the sequences of the inner vertices) and the goal is to reconstruct an approximation to the sequence of the root 
(the DNA of the ancestral father). For every e > 0, we present a deterministic algorithm which outputs an 
approximation of xi , . . . , Xk if pi + p<j < 0(1/ k 2 ^ 3 logn) and (1 — 2p s ) 2 > 0(d~ x log d). 

To our knowledge, this is the first rigorous trace reconstruction result on a tree in the presence of indels. 
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1 Introduction 



Trace reconstruction on a star. In the "trace reconstruction problem" (TRP) [LevOla, LevOlb, BKKM04, KM05, 
HMPW08, VS08], a random binary string X of length k generates an i.i.d. collection of traces Y±, . . . ,Y n that are 
identical to X except for random mutations which consist in indels, i.e., the deletion of an old site or the insertion 
of a new site between existing sites, and substitutions, i.e., the flipping of the state at an existing site 1 . (In keeping 
with biological terminology, we refer to the components or positions of a string as sites.) The goal is to reconstruct 
efficiently the original string with high probability from as few random traces as possible. 

An important motivation for this problem is the reconstruction of ancestral DNA sequences in computation 
biology [BKKM04, KM05]. One can think of X as a gene in an (extinct) ancestor species 0. Through speciation, 
the ancestor gives rise to a large number of descendants 1, . . . ,n and gene X evolves independently through 
the action of mutations into sequences Y\ , . . . , Y n respectively. Inferring the sequence X of an ancient gene from 
extant descendant copies Y±,. . . ,Y n is a standard problem in evolutionary biology [Tho04]. The inference of 
X typically requires the solution of an auxiliary problem, the multiple sequence alignment problem (which is an 
important problem in its own right in computational biology): site ti of sequence Yi and site tj of sequence Yj are 
said to be homologous (in this simplified TRP setting) if they descend from a common site t of X only through 
substitutions; in the multiple sequence alignment problem, we seek roughly to uncover the homology relation 
between Y\, . . . , Y n . Once homologous sites have been identified, the original sequence X can be estimated, for 
instance, by site-wise majority. 

The TRP as defined above is an idealized version of the ancestral sequence reconstruction problem in one 
important aspect. It ignores the actual phylogenetic relationship between species 1, . . . , n. A phylogeny is a 
(typically, binary) tree relating a group of species. The leaves of the tree correspond to extant species. Internal 
nodes can be thought of as extinct ancestors. In particular the root of the tree represents the most recent common 
ancestor of all species in the tree. Following paths from the root to the leaves, each bifurcation indicates a speciation 
event whereby two new species are created from a parent. An excellent introduction to phylogenetics is [SS03]. 

A standard assumption in computational phylogenetics is that genetic information evolves from the root to the 
leaves according to a Markov model on the tree. Hence, the stochastic model used in trace reconstruction can be 
seen as a special case where the phylogeny is star-shaped. (The substitution model used in trace reconstruction 
is known in biology as the Cavender-Farris-Neyman (CFN) [Cav78, Far73, Ney71] model.) It may seem that a 
star is a good first approximation for the evolution of DNA sequences. However extensive work on the so-called 
"reconstruction problem" in theoretical computer science and statistical physics has highlighted the importance of 
taking into account the full tree model in analyzing the reconstruction of ancestral sequences. 
The "reconstruction problem." In the "reconstruction problem" (RP), we have a single site which evolves through 
substitutions only from the root to the leaves of a tree. In the most basic setup which we will consider here, the 
tree is d-ary and each edge is an independent symmetric indel-free channel where the probability of a substitution 
is a constant p s > 0. The goal is to reconstruct the state at the root given the vector of states at the leaves. More 
generally, one can consider a sequence of length k at the root where each site evolves independently according to 
the Markov process above. Denote by n the number of leaves in the tree. The RP has attracted much attention in the 
theoretical computer science literature due to its deep connections to computational phylogenetics [Mos03, Mos04, 
DMR06, Roc08] and statistical physics [Mos98, EKPS00, MosOl, MP03, MSW04, JM04, BKMP05, BCMR06, 
GM07, BVVD07, Sly09a, Sly09c]. See e.g. [Roc07, Sly09b] for background. 

Unlike the star case, the RP on a tree exhibits an interesting thresholding effect: on the one hand, information is 
lost at an exponential rate along each path from the root; on the other hand, the number of paths grows exponentially 
with the number of levels. When the substitution probability is low, the latter "wins" and vice versa. This "phase 
transition" has been thoroughly analyzed in the theoretical computer science and mathematical physics literature — 
although much remains to be understood. More formally, we say that the RP is solvable when the correlation 
between the root and the leaves persists no matter how large the tree is. Note that unlike the TRP we do not require 

'One can also consider the case where X is arbitrary rather than random. We will not discuss this problem here. 
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high-probability reconstruction in this case as it is not information-theoretically achievable for d constant — simply 
consider the information lost on the first level below the root. Moreover the "number of traces" is irrelevant here as 
it is governed by the depth of the tree and the solvability notion implies nontrivial conelation for any depth. When 
the RP is unsolvable, the correlation decays to for large trees. The results of [BRZ95, EKPSOO, Iof96, BKMP05, 
MSW04, BCMR06] show that for the CFN model, if p s < p*, then the RP is solvable, where d(l - 2p*) 2 = 1. This 
is the so-called Kesten-Stigum bound [KS66]. If, on the other hand, p s > p* , then the RP is unsolvable. Moreover 
in this case, the correlation between the root state and any function of the character states at the leaves decays as 
n -^ 2 ( 1 ) The positive result above is obtained by taking a majority vote over the leaf states. 

Like the TRP, the RP is only an idealized version of the ancestral sequence reconstruction problem: it ignores 
the presence of indels. In other words, the RP assumes that the multiple sequence alignment problem has been 
solved perfectly. This is in fact a long-standing assumption in evolutionary biology where one typically prepro- 
cesses sequence data by running it through a multiple sequence alignment heuristic and then one only has to model 
the substitution process. This simplification has come under attack in the biology literature, where it has been 
argued that alignment procedures often create systematic biases that affect analysis [LG08, WSH08]. Much empir- 
ical work has been devoted to the proper joint estimation of alignments and phylogenies [TKF91, TKF92, Met03, 
MLH04, SR06, RE08, LG08, LRN+09]. 

Our results. We make progress in this recent new direction by analyzing the RP in the presence of indels — which 
we also refer to as the TRP on a tree (TRPT). We consider a d-axy tree where each edge is an independent channel 
with substitution probability p s , deletion probability p<±, and insertion probability p\ (see Section 1.1 for a precise 
statement of the model). The root sequence has length k and is assumed to be uniform in {0, l} h . As in the standard 
RP, we drop the requirement of high-probability reconstruction and seek instead a reconstructed sequence that has 
correlation with the true root sequence uniformly bounded in the depth. 

We give an efficient recursive procedure which solves the TRPT for p s > a small enough constant (strictly 
below, albeit close, to the Kesten-Stigum bound) and Pd,Pi = 0(k~ 2 / 3 log -1 n). As a by-product of our analysis 
we also obtain a partial global alignment of the sequences at the leaves. Our method provides a framework for 
separating the indel process from the substitution process by identifying well-preserved subsequences which then 
serve as markers for alignment and reconstruction (see Section 1.2 for a high-level description of our techniques). 
As far as we are aware, our results are the first rigorous results for this problem. 

Results on the RP have been used in previous work to advance the state of the art in rigorous phylogenetic tree 
reconstruction methods [Mos04, DMR06, MHR08, Roc08]. A central component in these methods is to solve the 
RP on a partially reconstructed phylogeny to obtain sequence information that is "close" to the evolutionary past; 
then this sequence information is used to obtain further structural information about the phylogeny. The whole 
phylogeny is built by alternating these steps. Our method sets up a framework for extending these techniques 
beyond substitution-only models. Partial results of this type will be given in the full version of the paper. 
Related work. Much work has been devoted to the trace reconstruction problem on a star [LevOla, LevOlb, 
BKKM04, KM05, HMPW08, VS08]. In particular, in [HMPW08], it was shown that, when there are only deletions, 
it is possible to tolerate a small constant deletion rate using poly(fc) traces. For a different range of parameters, 
Viswanathan and Swaminathan [VS08] showed that, under constant substitution probability and 0(1/ log k) indel 
probability, 0(log k) traces suffice. Both results assume that the root sequence X is uniformly random. 

The multiple sequence alignment problem as a combinatorial optimization problem (finding the best alignment 
under some pairwise scoring function) is known to be NP-hard [WJ94, Eli06]. Most heuristics used in practice, 
such as CLUSTAL [HS88], T-Coffee [NHH00], MAFFT [KMKM02], and MUSCLE [Edg04], use the idea of a 
guide tree, that is, they first construct a very rough phylogenetic tree from the data (using edit distance as a measure 
of evolutionary distance), and then recursively construct local alignments produced by "aligning alignments." Our 
work can be thought of as the first attempt to analyze rigorously this type of procedure. 

Finally, our work is tangentially related to the study of edit distance. Edit distance and pattern matching in 
random environments has been studied, e.g., by [NavOl, NBYST, AK08]. 
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1.1 Definitions 



We now define our basic model of sequence evolution. 

Definition 1.1 (Model of sequence evolution) Let be the d-ary tree with H levels and n = d H leaves. For 

simplicity, we assume throughout that d is odd. We consider the following model of evolution on Tjp. The sequence 

at the root ofT^ has length k and is drawn uniformly at random over {0, l} k . Along each edge of the tree, each 
site (or position) undergoes the following mutations independently of the other sites: 

• Substitution. The site state is flipped with probability p s > 0. 

• Deletion. The site is deleted with probability p^ > 0. 

• Insertion. A new site is created to the right of the current site with probability p\ > 0. The state of this new 
site is uniform {0, 1}. 

These operations occur independently of each other. The last two are called indels. We let p\d = p\ + Pd and 
6 S = 1 — 2p s . The parameters p s ,Pd,Pi may depend on k and n, where n is the number of leaves. 

Remark 1.2 For convenience, our model of insertion is intentionally simplistic. In the biology literature, related 
continuous-time Markov models are instead used for this kind of process [TKF91, TKF92, Met03, MLH04, RE08, 
DR09 ]. It should be possible to extend our results to such generalizations by proper modifications to the algorithm. 

1.2 Results 

Statement of results. Our main result is the following. Denote by X = x\, . . . , x^ a binary uniform sequence of 
length k. Run the evolutionary process on with root sequence X and let Y\ , . . . , Y n be the sequences obtained 
at the leaves, where Yi = y\, . . . , . 

Theorem 1.3 (Main result) For all x > and (5 = 0(d~ 1 log d), there is <3?, 3>', <J>" > such that the following 
holds for d large enough. There is a polynomial-time algorithm A with access to Yi, . . . ,Y n such that for all 

(l-2p s ) 2 >^^, •p i+ p d <_Jj , <£" log 3 n < k < poly(n), 

the algorithm A outputs a binary sequence X which satisfies the following with probability at least 1 — x : 

1. X = x\, . . . , Xk has length k. 

2. For allj = l,..., k, P[% = xj] > 1 - (3. 

Remark 1.4 Notice that we assume that the (leaf-labelled) tree and and the sequence length of the root are known. 
The requirement that the sequence length is known is not crucial. We adopt it for simplicity in the presentation. 

Remark 1.5 In fact, we prove a stronger result which allows x = o(l) and shows that the "agreement" between 
X and X "dominates" an i.i.d. sequence. See Lemma B.l and Section 5.2. 

Proof sketch. We give a brief proof sketch. As discussed previously, in the presence of indels the reconstruction 
of ancestral sequences requires the solution of the multiple sequence alignment problem. However, in addition to 
being computationally intractable, global alignment through the optimization of a pairwise scoring function may 
create biases and correlations that are hard to quantify. Therefore, we require a more probabilistic approach. From 
a purely information-theoretic point of view the pairwise alignment of sequences that are far apart in the tree is 
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difficult. A natural solution to this problem is instead to perform local alignments and ancestral reconstructions, 
and recurse our way up the tree. 

This recursive approach raises its own set of issues. Consider a parent node and its d children. It may be easy to 
perform a local alignment of the children's sequences and derive a good approximation to the parent sequence (for 
example, through site-wise majority). Note however that, to allow a recursion of this procedure all the way to the 
root, we have to provide strong guarantees about the probabilistic behavior of our local ancestral reconstruction. 
As is the case for global alignment, a careless alignment procedure creates biases and correlations that are hard to 
control. For instance, it is tempting to treat misaligned sites as independent unbiased noise but this idea presents 
difficulties: 

Consider a site j of the parent sequence and suppose that for this site we have succeeded in aligning 
all but two of the children, say 1 and 2. Let x\ denote the site in the z'th child which was used to 
estimate the j'th site. By the independence assumption on the root sequence and the inserted sites, xj x 
and x? 2 are uniform and independent of (x* .)f =3 . However, x 1 , and xj 2 may originate from the same 
neighboring site of the parent sequence and therefore are themselves correlated. 

Quantifying the effect of this type of correlation appears to be nontrivial. 

Instead, we use an adversarial approach to local ancestral reconstruction. That is, we treat the misaligned 
sites as being controlled by an adversary who seeks to flip the reconstructed value. This comes at a cost: it 
produces an asymmetry in our ancestral reconstruction. Although the RP is well-studied in the symmetric noise 
case, much remains to be understood in the asymmetric case. In particular, obtaining tight results in terms of 
substitution probability here may not be possible as the critical threshold of the RP may be hard to identify. We do 
however provide a tailored analysis of the particular instance of the RP by recursive majority obtained through this 
adversarial approach and we obtain results that are close to the known threshold for the symmetric case. Unlike 
the standard RP, the reconstruction error is not i.i.d. but we show instead that it "dominates" an i.i.d. noise. (See 
Section 4.2 for a definition.) This turns out to be enough for a well-controlled recursion. We first define a local 
alignment procedure which has a fair success probability (independent of ri). However, applying this alignment 
procedure multiple times in the tree is bound to fail sometimes. We therefore prove that the local reconstruction 
procedure is somewhat robust in the sense that even if one of the d inputs to the reconstruction procedure is faulty, 
it still has a good probability of success. 

As for our local alignment procedure, we adopt an anchor approach. Anchors were also used by [KM05, 
HMPW08] — although in a quite different way. We imagine a partition of every node's sequence into islands of 
length Oik 1 / 3 ). (The precise choice of the island length comes from a trade-off between the length and the number 
of islands in bounding the "bad" events below — see the proof of Lemma 3.3.) At the beginning of each island 
we have an anchor of length O (log re). Through this partition of the sequences in islands and anchors we aim to 
guarantee the following. Given a specific father node v, with fair probability 1) all the anchors in the children 
nodes are indel-free; and 2) for all parent islands, almost all of the corresponding children islands have no indel at 
all and, moreover, at most one child island may have a single indel. The "bad" children islands — those that do not 
satisfy these properties — are treated as controlled by an adversary. We show that Conditions 1) and 2) are sufficient 
to guarantee that: the anchors of all islands can be aligned with high probability and single indel events between 
anchors can be identified. This allows a local alignment of all islands with at most one "bad" child per island and is 
enough to perform a successful adversarial recursive majority vote as described above. The bound on the maximum 
indel probability sustained by our reconstruction algorithm comes from satisfying Conditions 1) and 2) above. 
Notation. For a sequence X = x\, . . . , x^, we let X[i : j] = aij, . . . , Xj. We use the expression "with high 
probability (w.h.p.)" to mean "with probability at least 1 — l/poly(n)" where the polynomial in n can be made 
of arbitrarily high degree (by choosing the appropriate constants large enough). We denote by Bin(ra, p) a random 
variable with binomial distribution of parameters n,p. For two random variables X, Y we denote by X ~ Y the 
equality in distribution. 

Organization. The rest of the paper is organized as follows. We describe the algorithm in Section 2. The proof of 
our main result is divided into two sections. In Section 3, we prove a series of high-probability claims about the 
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evolutionary process. Then, conditioning on these claims, we provide a deterministic analysis of the correctness of 
the algorithm in Section 5.2. All proofs are in the Appendix. 

2 Description of the Algorithm 

In this section we describe our algorithm for TRPT. Our algorithm is recursive, proceeding from the leaves of the 
tree to the root. We describe the recursive step applied to a non-leaf node of the tree. 

Recursive Setup — Our Goal. For our discussion in this section, let us consider a non-leaf node v with d children, 
denoted m for i £ [d]. For notational convenience, we drop the index u and denote its children by 1, . . . , d. Our 
goal for the recursive step of the algorithm is to reconstruct the sequence at the node v given the sequences of the 
children. Denote the sites of the father by Xo = Xi,..., x^ o , and the sites of the i'th child by Xi = x\ , . . . , x\.. 
During the reconstruction process, we do not have access to the children's sequences, but rather to reconstructed 
sequences denoted by Xj — X .. . . . X * . 

Let us consider the following partition of the sequence of v into subsequences, called islands. Of course our 
algorithm doesn't have access to the sequence at v during the recursive step of the algorithm. We define the partition 
as a means to describe our algorithm: The sites of v are partitioned into islands of length £ = /c 1 / 3 (except for the 
last one which is possibly shorter). Denote by iVo = \ko/£\ the number of islands in v. Each island starts with 
an anchor of a bits. That is, the islands are the bitstrings Xq\1 : £], Xq[£ + 1 : 2£], . . . and the anchors are the 
bitstrings X [l : a], X [£ + 1 : £ + a], . . .. 

Our algorithm tries to identify for each island Xo[(i — 1)1 + 1 : i£] the substrings of each of the d children 
that correspond to this island (i.e., contain the sites of the island), called "child islands." We do so iteratively for 
i = 1 . . . iVo- We use the islands that did not have indels for sequence reconstruction, using the substitution-only 
model. Some islands will have indels however. This leads to two "modes of failure": one invalidates the entire 
(parent) node, and the other invalidates only an island of a child. More specifically, a node becomes invalidated 
(i.e., useless) when indels are not evenly distributed, that is: when an indel occured in an anchor, or two (or more) 
indels occured in a specific island over all d children. This is a rare event. Barring this event, we expect that 
each island suffers only at most one indel over all children. The island (of a child) that has exactly one indel is 
invalidated (second mode of failure), and is thus deemed useless for reconstruction puiposes. As long as the parent 
node is not invalidated, each island will have at least d — 2 non-invalidated children islands (one additional island 
is potentially lost to a child node that may have been invalidated at an earlier stage). 

Even when the algorithm identifies that a child island has an indel somewhere, the island is not ignored. The 
algorithm still needs to compute the length of the island in order to know the start of the next island in this child. 
For this purpose, we use the anchor of the next island and match it to the corresponding anchors of the other 
(non-invalidated) child islands. In fact the same procedure lets us detect which of the child islands are invalidated. 

More formally, we define d functions /j : {1, . . . , ko} — ► {l,...,fcj}U{f}, where /j takes the sites of v to 
the corresponding sites of the i'th child or to the special symbol f if the site was deleted. Note that for each i, /j 
is monotone, when ignoring sites which are mapped to f. For t = £r, let Sj(r) = fi(t + 1) — (t + 1) denote the 
displacement of the site corresponding to the (t + l) st site of the parent, in the i th child. By convention, we take 
Si(0) = 0. If there is no indel between t = £r and t' = £r' then Si(r) = Sj(r'). Note that, in the specific case of 
one indel operation in the island, we have that |sj(r) — Sj(r')| = 1. 

Algorithm. Our algorithm estimates the values of Sj(r) and uses these estimates to match the stalling positions 
of the islands in the children. The full algorithm is given in Figure 1 in the Appendix. We use the following 
additional notation. For x € {0, 1} we let (x) = 2x — 1. Then, for two {0, l} m -sequences Y = yi, ■ ■ ■ ,y m and 
Z = z\, . . . , z m , we define their (empirical) correlation as 

1 m 

CoTv(Y,Z) = -^2(y j )(z j ). 
Note that y t— > (y) maps 1 to 1 and to —1. One can think of Corr(Y, Z) as a form of normalized centered 
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Hamming distance between Y and Z. In particular, a large value of Corr(Y, Z) implies that Y and Z tend to agree. 
We will use the following threshold (which will be justified in Section 5.1) 

7 = ((l-5)(l-2 Ps ) 2 -4/3), 

where 5 is chosen so that 

{l-6)(l-2 Ps ) 2 -8f3> 8 + 8(3, 

where again j3 = 0(d~ l log d). 

3 Analyzing the Indel Process 

We define a > Clogn and a < e/d < 1, for constants C, e to be determined later. We require a < fc 1 / 3 < 
poly(n). 2 We assume that the indel probability per site satisfies 



1 4dk 2 / 3 a \k 2 / 3 lognJ 
Throughout, we denote the tree by T = (V, E). 

3.1 Bound on the Sequence Length 

As the indel probability is defined per site, longer sequences suffer more indel operations than shorter ones. We 
begin by bounding the effect of this process. We show that with high probability the lengths of all sequences are 
roughly equal. 

Lemma 3.1 (Bound on sequence length) For all ( > (small), there exists C > (large) so that for all u in V, 

we have 

k v G [k, k] = [(1 — (1 + C)^L 
with high probability given k > C log 3 n. We denote this event by C 

3.2 Existence of a Dense Stable Subtree 

In this section, we show that with probability close to 1 there exists a dense subtree of T with a "good indel 
structure," as defined below. Our algorithm will try to identify this subtree and perform reconstruction on it, as 
described in Section 4. 

Indel structure of a node. Recall that £ = fc 1 / 3 . 

Definition 3.2 (Indel structure) For a node (parent) v, we say that v is radioactive if one of the following events 
happen: 

1. Event B\: Node v has a child u such that when evolving from v to u an indel operation occurred in at least 
one of the sites which are located in an anchor. 

2. Event B%: There is an island I and two children u, u', such that an indel occurred in I in the transition from 
v to u and in the transition from v to v! . 

3. Event B3: There is an island I and a child u, such that two indel operations (or more) happened in I in the 
transition from v to u. 

2 A variant of the algorithm where the anchors have length 0(log k) also works when k 2> n. 



6 



Otherwise the node v is stable. By definition, the leaves ofT are stable. A subtree ofT is stable if all of its nodes 
are stable. 

Lemma 3.3 (Bound on radioactivity) For all a > 0, there exists a choice of Q > small enough in Lemma 3.1 
such that conditioning on the event C occuring: any vertex v is radioactive with probability at most a. 

As a corollary we obtain the following. 

Lemma 3.4 (Existence of a dense stable subtree) For all x > 0, there is of ( > small enough in Lemma 3. 1 
such that, conditioning on the event C occuring, with probability at least 1 — x> the root of T is the father of a 
(d — l)-ary stable subtree ofT. We denote this event by S. 

4 A Stylized Reconstruction Process 

In this subsection, we lay out the basic lemmas that we need to analyze our ancestral reconstruction method. We 
do this by way of describing a hypothetical sequence reconstruction process performed on the stable tree defined 
by the indel process (see Lemma 3.4). We analyze this reconstruction process (assuming that the radioactive 
nodes and the islands with indels are controlled by an adversary) and show in Lemma 4.5 that the process gives 
strong reconstruction guarantees. Then we argue in Section 5 that our algorithm performs at least as well as the 
reconstruction process against the adversary described in this section. Throughout this section we suppose that a 
stable tree exists and is given to us, together with the "orbit" of every site of the sequence at the root of the tree (see 
function F below). However, we are given no information about the substitution process. 

Let v G V and assume v is the root of a (d — l)-ary stable subtree T* = (V*, E*) of T. (We make the stable 
subtree below v into a (d— l)-ary tree by potentially removing arbitrary nodes from it, at random.) Let u G V*. For 
each island I in u, at most one child v! of u in T* contains an indel in which case it contains exactly one indel. We 
say that such an / is a corrupted island of v! . The basic intuition behind our analysis is that, provided the alignment 
on T* is performed correctly (which we defer to Section 5.2), the ancestral reconstruction step of our algorithm is 
a recursive majority procedure against an adversary which controls the corrupted islands and the radioactive nodes 
(as well as all their descendants). Below we analyze this adversarial process. 

Recursive majority. We begin with a formal definition of recursive majority. Let Maj : {0, 1, — > {0, 1} be the 
function that returns the majority value over non-ft values, and flips an unbiased coin in case of a tie (including the 
all-(J vector). Let no = d H ° be the number of leaves in T below v. Consider the following recursive function of 

z = (z 1 ,z 2 ,..., z no ) G {0, 1, Jt}: Maj (zi) = Z\, and 

Maj J (zi, ...,z dJ ) = Maj(Maj : '" 1 (2ii, . . . , z d u-i)), Maf -1 ^.^--!)^, . . . , z dJ )), 

for all j = 1, ... , Hq. Then, M&i H °(z) is the d-wise recursive majority of z. 

Let Xq = Xj 1 , • • • , x\ Q be the sequence at v. For u G V* and t = 1, . . . , ko, we denote by F u (t) the position of 
site x® in u or f if the site has been deleted on the path to u. We say that C u j holds if F u (t) is in a corrupted island 
of u. Let Path(n, v) be the set of nodes on the path between u and v. 

Definition 4.1 (Gateway node) A node u is a gateway for site t if: 

L F u (t) / f; and 

2. For all u' G Path(u, v) — {v}, C u i ± does not hold. 

We let = (V t **,E**) be the subtree of T* containing all gateway nodes for t. By construction, is at least 
(d — 2)-ary and for convenience we remove arbitrary nodes, at random, to make it exactly (d — 2)-ary. Notice that, 
for t,t' G [1 : ko], the subtrees T t ** and T t ** are random and correlated. However, they are independent of the 
substitution process. 
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We will show in Section 5.2 that the reconstructed sequence produced by our method at v "dominates" (see 
below) the following reconstruction process. Let L v = u\ , . . . , u no be the leaves below v ordered according to a 
planar realization of the subtree below v. Denote by Xj = x\, . . . , x\. the sequence at For t = 1, . . . , ko, let 
L\* be the leaves of T 4 **. We define the following auxiliary sequences: for m € L v , we let 3j = £f , . . . , £L where 
for t = 1, . . . , &o 

1 — x" o.w. 

In words, £| is the descendant of if u\ is a gateway to t and is the opposite of the value x® otherwise. Because 
of the monotonicity of recursive majority, the latter choice is in some sense the "worst adversary" (ignoring cor- 
relations between sites — we will come back to this point later). We then define a reconstructed sequence at v as 
3o = ii, ■ ■ ■ , Ck where for £ = 1, . . . , fco 

We now analyze the accuracy of this (hypothetical) estimator — which we refer to as the adversarial reconstruction 
of Xq. We show in Section 5.2 that our actual estimator is at least as good as 3o w.h.p. 

4.1 Recursive Majority Against an Adversary 

To analyze the performance of the adversarial reconstuction So, we consider the following stylized process. 
Definition 4.2 (Recursive Majority Against an Adversary) We consider the following process: 

1. Run the evolutionary process on T^ o ^ at one position only starting with root state without indels, that is, 
taking p 1( \ = 0. 

2. Then complete Tjj o 2 ^ into and associate to each additional node the state 1. 

3. Let Rfl o be the random variable in {0,1} obtained by running recursive majority on the leaf states obtained 
above. 

We call this process the recursive majority against an adversary on . We show the following. 

Lemma 4.3 (Accuracy of recursive majority) For all > 0, there exists a constant C" > such that taking 

and d large enough, then the probability that the recursive majority against an adversary on Tg correctly recon- 
structs root state is at least 1 — (3 uniformly in Hq. In comparison, note that the Kesten-Stigum bound for binary 
symmetric channels on d-ary trees is 9 2 > d" 1 [KS66, Hig77]. 

As a corollary of Lemma 4.3, we have the following. 

Definition 4.4 (Bernoulli sequence) For q > and m G N, the (q, m)-Bernoulli sequence is the product dis- 
tribution on {0, l} m such that each position is 1 independently with probability 1 — q. We denote by B„ m the 
corresponding random variable. 

Lemma 4.5 (Subsequence reconstruction) Assume v is the root of a (d — \)-ary stable subtree. For all (3 > 0, 

choosing C" > as in Lemma 4.3 is such that the following holds for d large enough. For t,m G {1, . . . , ko}, let 
A = (Ai, . . . , A m ) be the agreement vector between the Eo[t + 1 : t + m] and Xo[t + 1 : t + rn], that is, Aj = 1 if 
recursive majority correctly reconstructs position i. Then there is < [3 1 < (3 such that A ~ Bp m . (Here, may 
depend on Hq but (3 does not.) 
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4.2 Stochastic Domination and Correlation 



In our discussion so far we have assumed that a stable tree exists and is given to us, together with the the function 
F. This allowed us to define the stylized recursive majority process against an adversary (Definition 4.2), for which 
we established strong reconstruction guarantees (Lemmas 4.3 and 4.5). In reality, we have no access to the stable 
tree. We are going to construct it recursively from the leaves towards the root. At the same time we will align 
sequences, discover corrupted islands, and reconstruct sequences of internal nodes. The stylized recursive majority 
process will be used to provide a "lower bound" on the actual reconstruction process. The notion of "lower bound" 
that is of interest to us is captured by stochastic domination, which we proceed to define formally. 

Definition 4.6 (Stochastic domination) Let X,Y be two random variables in {0, l} m . We say that Y stochasti- 
cally dominates X, denoted X < Y, if there is a joint random variable (X,Y) such that the marginals satisfy 
X ~ X andY ~ Y and moreover F[X < Y] = 1. 

Note that in the definition above X and Y may (typically) live in different probability spaces. Then, the joint 
variable (X, Y) is a coupled version of X and Y. In our case, X is the adversarial recursive process whereas Y 
is the actual reconstruction performed by the algorithm. We now show how to use this property for correlation 
estimation. 

Correlation. The analysis of the previous section guarantees that the sequences output by the adversarial recon- 
struction process are well correlated with the true sequences. But if we are only going to use the adversarial 
process as a lower bound for the true reconstruction process, it is important to establish that stochastic domina- 
tion preserves correlation. In preparing the ground for such a claim let us establish an important property of the 
adversarial process. Let T u and T v be the two disjoint copies of rooted at the nodes u and v respectively, 
and let X = xi, X2, ■ ■ ■ , x m G {0, l} m and Y = y±, y2, ■ ■ ■ , y m G {0, l} m be sequences at the nodes u and 



v. 



Assume that u and v are the roots of (d — l)-ary stable subtrees. Let X' 



ry' ry' 



G {0, l} m and 



Y' = y[, y' 2 i ■ ■ • > y'm ^ {0> l} m De the reconstructions of X and Y obtained by the adversarial reconstruction 
process. Let A = Ai, . . . , A m and = 9±, . . . , 6 m be the resulting agreement vectors. We show the following: 



Lemma 4.7 (Concentration of bias) Let (3', (3 be as in Lemma 4.5. Then, with probability at least 1 
the following are satisfied 



1 E(A 1 )(^)-(l-2/? / ) 2 



^ m 



i=l 



i=l 



a, 



We use the previous lemma to argue that stochastic domination does not affect our correlation computations. 

Lemma 4.8 (Correlation bound) Let X, Y G {0, l} m be random strings defined on the same probability space 
as X' and Y'. Denote by Z (resp. W) the agreement vectors of X (resp. Y) with X (resp. Y). Assume that A < Z 
and G < W with probability 1, where A and are the agreement vectors of X' and Y' with X and Y as explained 
above. Then, 

^ m 

| Corr(A, Y) - Cott(X,Y)\ < 1 - - V((Ai)(^ 



m 



i=l 



1 



<Ai>= 



1/ 



k)- 



with probability 1. Furthermore, conditioned on the conclusions of Lemma 4.7, we have, with probability 1: 



Corr(A, Y) - Corr(X, Y)| < 
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5 Analyzing the True Reconstruction Process 



We provide the proof of Theorem 1.3. In Section 5.1, we show that, if a stable subtree exists, the adversarial 
reconstructions of aligned anchors exhibit strong correlation signal, while misaligned anchors exhibit weak signal. 
This holds true for sequences that stochastically dominate the adversarial reconstructions. We use this property to 
complete the analysis of our reconstruction method in Section 5.2. 



5.1 Anchor Alignment 



Consider a parent v that is stable. Let i, j be two children with sequences Xi = x\, . . . , x\_ and Xj = x\, . . . , x 3 k 
Let t = tr and consider the following subsequences (of length a) at i and j 

8% = x i [t + Si(r) + 1 : t + Si(r) + a], and satf = x j [t + sj(r) + 1 : t + Sj(r) + a]. 

These are related (but not identical) to the definition of anchors in the algorithm of Section 2. In particular, note 
that by definition and are always aligned, in the sense that they correspond to the same subsequence of v. 
Consider also the following subsequences 

@ 3 r = x j [t + Sj(r) : t + Sj(r) + a - 1] and = x j [t + Sj(r) + 2 : t + Sj(r) + a + 1]. 

These are the one-site shifted subsequences for j. The following lemma bounds the correlation between these 
strings. More precisely, we show that is always significantly more correlated to its aligned brother gf$ than 
to the misaligned ones *2)i and JP£ . This follows from the fact that the misaligned subsequences are sitewise 
independent. 

Lemma 5.1 (Anchor correlations) For all 5 > such that (1 — 5) (1 — 2p s ) 2 — 8/3>5 + 8/3, there is C > large 
enough so that with a = Clog re, the following hold: 



1. Aligned anchors. 



Corr^;, s4) > (1 - 6)(1 - 2 Ps ) 2 > 1 - exp (-0(a)) = 1 - l/poly(n) 



2. Misaligned anchors. P 



Corr(<<,^) < 5 



> 1 — exp ( — 0(a)) = 1 — l/poly(re), and similarly for ,_S?r ■ 



We denote by Aij r the above events and their symmetric counterparts under i <-» j. 

Lemma 5.2 (Anchor correlations: Reconstructed version) Let Xi = {x\) k L L x and Xj = (xi)^ dominate the 
adversarial reconstructions X- and Xj of Xi and Xj, as defined in Lemma 4.8. Let ^ = x l [t + Sj(r) + 1 : 
t + Sj(r) + a) and similarly for all other possibilities srf «-> J? and/or i <-> j. Denote by Bij r the event that the 
conclusions of Lemma 4.7 hold for X^ and X'- over all pairs of intervals involving [t + Sj(r) : t + Sj(r) + a — 1], 
[t + Si(r) + 1 :t + Si(r) + a], and [t + Sj(r) + 2 : t + s»(r) + a + 1], with i <-> j <xs necessary. Then, conditioned 
on Bij t1 . and Aij >r we have 

Corr(<,< i ) > (1 " <*)(! " 2p s ) 2 - 8/3, 
Corr(^;, ^') < 5 + 8/3, Corr(^, ^') < 5 + 8/3, 
we/Z their symmetric counterparts under i <-> j. 



5.2 Proof of Correctness 

We show that our recursive procedure reconstructs the desired sequence at the root of the tree whenever a collection 
of good events occurs. Recall the definitions of the events C, S, Bi,j jT , Ai.j jT from Lemmas 3.1, 3.4, 5.1 and 5.2. 3 

3 Event £ guarantees that there is no big variance in the nodes' sequence lengths; event S guarantees that a stable (d — l)-ary subtree 
exists; the events Bi t j. r guarantee that the adversarial reconstruction process is successful, also in preserving correlations between sequences 
of nodes; and the events Ai.j.r guarantee that aligned anchors (across sequences of a node's children) exhibit strong correlation signal, while 
misaligned anchors give weak correlation signal. 
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Conditioning on C and S, denote by T* = (V*,E*) the stable (d - l)-ary subtree of T. Then, for all v G V* , 
all pairs of children i,j of v in T*, and all r = 1, . . . , we condition on the events Bi^ r and Ai,j ;r - Note that 
having conditioned on C there is only a polynomial number of such events, since all sequence lengths are bounded 
by k. (If r ■ t is larger than a node's sequence length we assume that the corresponding events are vacuously 
satisfied.) Finally recall that, conditioning on C, the event S occurs with probability 1 — x an d all other events 
occur with high probability. We denote the collection of events by £. 

Conditioning on £, the proof of correctness of the algorithm follows from a bottom-up induction. The gist of 
the argument is the following. Suppose that at a recursive step of the algorithm we have reconstructed sequences 
for all children of a node v, which are strongly correlated with the true sequences (in the sense of dominating the 
corresponding adversarial reconstructions). Having conditioned on the events Ai,j, r and Bij jT , it follows then that 
the correct alignments of anchors exhibit strong correlation signal while the incorrect alignments weak correlation 
signal. Hence, our conelation tests between anchors discover the corrupted islands and do the anchor alignments 
correctly (at least for all nodes lying inside the stable tree). Hence the shift functions Sj's are correctly inferred, 
and the reconstruction of v's sequence can be shown to dominate the corresponding adversarial reconstruction. The 
complete proof details are given in App endix D. 
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A Algorithm 



1. Input. Children sequences x 1 , . . . ,x d . 

2. Initialization. Set sj(0) := 0, Vi, £ = fc 1/3 , r = 1, and t = I 

3. Main loop. While x l [t + Sj(r — 1) + 1 : i + Sj(r — 1) + a] is non-empty for all i, 

(a) Current position. Set t = £r. 

(b) Anchor definition. For each i, set A\ = x l [t + §i(r — 1) + 1 : t + §i(r — 1) + a]. 
We say that A\ is the r'th anchor of the z'th child. (If the remaining sequences are not 
long enough to produce an anchor of length a, we repeat the previous step with the full 
remaining sequences.) 

(c) Alignment. For each anchor, we define the set of anchors which agree with it. Formally, 

G; = 0e[4Corr(l*,^)> 7 }. 

(d) Update. We define the set of aligned children G r = {i : \G l r \ > d — 2}. 

i. Aligned anchors. For each i G G r , set Sj(r) = ij(r — 1). 

ii. Misaligned anchors. For each i £ G r define two strings D l r = x l [t + 5j(r — 1) : 
t + Si(r - 1) + a - 1] and T* = x i [t + §i(r - 1) + 2 : t + s,(r - 1) + a + 1]. If 

10 G [41 " W : Corr(5*,J4?) > 7 }| > d- 2, 

set Sj(r) = Sj(r — 1) — 1. If 

\{je[d]-{i} : Corr(? r ,^)> 7 }|>d-2, 

set §i(r) =§i(r — l) + 1. 

(e) Ancestral sequence. Compute ... x^ 1 by performing a sitewise majority on the 
children in G r . (If the remaining children sequences are too short to produce a full 
island, we use whatever is left which should all have equal length by our proof.) 

(f) Increment. Set r := r + 1. 

4. Output. Output x° and set ko to its length. 

Figure 1 : This is the basic recursive step of our reconstruction algorithm. It takes as input the d inferred sequences 
of the children x l , . . . , x d and computes a sequence for the parent x°. If any of the steps above cannot be accom- 
plished, we abort the reconstruction of the parent and declare it radioactive. 

B Further Lemmas 

For a going to 0, we have more precisely: 
Lemma B.l (Limit a — > 0) Condition on C. Let 

a 
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for h(n) = to(l). Then, for n large enough, the root is the father of a (d — l)-ary stable subtree with probability at 
least 

1-X= i — 

y/Hn) 

Proof of Lemma B.l: Plugging a = l/h{n) and v = 1 — h{n) into the recursion derived in the proof of 
Lemma 3.4, we get 

(l-a)g(v) = (l-J—\ll-—£ = + o( T ^) 
V h(n)J \ ^fh{n) \h(n)J 

> l-i/v^H, 



C Proofs 

Proof of Lemma 3.1: We prove the upper bound by assuming there is no deletion. The lower bound can be proved 
similarly. The proof goes by induction. Let v be a node at graph distance i from the root. We show that there is 
C" > independent of i such that 

k v <k + i^C'klogn. 
Since the depth of T is 0(log n), this implies the main claim as long as 



\J C"k log n log n < (k, 

which follows from our assumption for C > large enough. 

The base case of the induction is satisfied trivially. Assume the induction claim holds for v, the parent of u. It 
suffices to show that the number of new insertions is at most \J C"k log n. By our induction hypothesis, the number 
of insertions is bounded above by a binomial Z with parameters k + (i — l)yJC"k log n < (1 + Qk and w.h.p. 
By Hoeffding's inequality, taking 

IC" logn 



(i + C)*' 



we have 



F[Z > (1 + ()kp id + (1 + 0*7/] < exp(-2((l + C)A:r ? ) 2 /[(l + C)fc]) 

= Vpoly(ra). 

By our assumption on p\&, we have 

( ak l l^\ 

(1 + C)k m = O 



log n 



so that choosing C" large enough gives 



(1 + C)k Pid + (1 + ()k V < y/C"k\og 



n. 
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This proves the claim. ■ 



Proof of Lemma 3.3: According to Lemma 3.1, the length of the sequence at v is in [k, k] w.h.p. We denote that 
event by C v . We bound the probability of events 81,82, B3 separately. 

Let N = kjl = (1 + C)& 2 ^ 3 - Conditioned on C v , there are at most N anchors, each of length a. By a union 
bound, the probability that at least one of the sites in the anchors has an indel operation in any child is upper 
bounded by 

P[Bi] = P[£i| +P[£i|£S]P[££] 
< Nadpid + l/poly(n) 
aadN 



Ak 2 / 3 ad 



= k 2/s ■ 4 + Vpoiy(n) 

< a(l/3 - l/poly(n)), 

where we choose ( small enough. The quantity we want to estimate is in fact P[B± | C] (which is not the same as 
conditioning on C v only). But notice that 

P[Bi] = F[Bi I C]F[C] + P[Bi I £ C ]P[£ C ] > P[B X | £]P[£], 

which implies 

mrR l n ^ Q(V3 - l/poly(n)) 

P[Bl 1 C] ~ 1 - l/poly(n) < a/3 ' 

(This argument shows that it suffices to condition on C v . We apply the same trick below.) 

To bound the probability of the second event, consider an island I and a son u. The probability that there is an 
indel when evolving from v to u is at most 

Ad* — , ,9/, _ 



4k 2 / 3 ad Ak l /z a d' 
Thus, the probability that more than one child of v experiences an indel in I is at most 

d 



Efd\ ( a y ^ d l / a y 
I J UfcV3 a J " f-f if UfcV3 a J 



i=2 v x i=2 

J 



< 



El / a 
i! UfcVSa 

i=2 



a 

< e 



4A; 1 /3 a 
ea 2 
16A; 2 / 3 a 2 ' 

where we used that the expression in parenthesis on the second line is < 1. Taking a union bound over all islands, 
the probability that at least two children experience an indel in the same island is at most 



F[B 2 1 C] < N 



ea 2 



16/t 2 / 3 o 2 
1 + C)ea 2 



16a 2 



a 

< 3' 
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where we used that a < 1. 

For the third event, consider again an island / and a child u. The probability for at least two indel operations in 
I when evolving from v to u is at most 



21 

£ 

i=2 



2( 



a 



4adk 2 /z 



< 



< 



1 / 2£a V 



i=2 
2t 



5^ i! 



a 



i=2 

< e 
< 



a 



2adfcV3 

2 



2adk 1 / 3 
ea 2 



Aa 2 d 2 k 2 / 3 ' 

(We use 21 to account for insertions and deletions.) Taking a union bound over all islands and children, the 
probability that there are two indel operations in the same child in the same island is bounded by 

2 

F[B 3 \C] < dN 



< 



4a 2 d 2 k 2 / 3 
(1 + C)e« 2 



Aa 2 d 



< a/3. 

Taking a union bound over the three ways in which a site can become radioactive proves the lemma. ■ 

Proof of Lemma 3.4: We follow a proof of [MosOl]. Let v be a node at distance r from the leaves. We let v r be 
the probability that v is the root of a (d — l)-ary stable subtree. Let 

d-l. 



Then, from Lemma 3.3, 
Note that 



g [v) = v d + dv d -\\ - v). 



v r > (1 - a)g(u r -i). 



g'(u) = d(d-2)u d - 2 (l-u). 
In particular, g is monotone, g(l) = 1, and g'(l) = 0. Hence, there is 1 — \ < v * < 1 su °h that 

g[y )>v . 

Then, taking 



1 - a > v*/g(v*), 



we have 



v r > (1 - a)g(v r _i) > 



g(y*) 



g{v r -\) > > i - x, 



by the induction hypothesis that v r ~i > v*- Note in particular that vq = 1 > v* . ■ 

Proof of Lemma 4.3: Recall that we assume the root state is and all adversarial nodes are 1. Because of the 
bias towards 1, we cannot apply standard results about recursive majority for symmetric channels [Mos98, Mos04]. 
Instead, we perform a tailored analysis of this particular channel. 
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We take asymptotics as d — > +00 and we show that the probability of reconstruction can be taken to be 



I -" = 1 -3- 



for C" large enough. Let v be the root of Tg . We denote by Z v the number of non-adversarial children of v in state 
and by Z' v the number of nodes among them that return upon applying recursive majority to their respective 
subtree. Let q Q HQ be the probability of incorrect reconstruction at v (given that the state at v is 0). Then 



> 



> 



Z'> 



d+l 



d-2 

E 

i=0 



Z 1 ,, > — - — I Z v 



nz v 



a) 



where we simply ignored the contribution of the children who flipped to 1. 

We prove q° Ho < 1/d by induction on the height. Let u be a non-adversarial node in at height h from the 
leaves to which we associate as above the variables Z u , Z' u and the quantity gP. Note that q$ = 0. We assume the 
induction hypothesis holds for h — 1. Note that conditioned on the state at u being Z u is Bin(d — 2, (1 — p s )) 
where 

l + #s 



1 - Ps 



+ e 



logd 
d 



as d — > +00. Similarly, given Z M = i, the variable is Bin(z, 1 — q^-i)- in particular, the quantity 

z> > d + l \z -i 

is monotone in i We use Chernoff 's bound on Z' u to truncate the lower bound (1). Indeed, let 

d 



M=(l-ft)(d-2) = - + T(d) ) 



with 



and 



T(d) = e(^dlogd), 



x d T(d) 
M(l-^) = 2 + 4 i 



where in particular 



Then, we have 



log d 



F[Z U < - t?)] < exp (- W 2 /2) = rf-nW, 
for C" large enough. Applying to (1) leads to the lower bound 



l- q ° h >(l-d-^] 



Bin 



d m \ d+i 

2 + 2 , 1 I > 2 



By the induction hypothesis, q9 1 < By applying Chernoff 's bound again we get 



. d T d 
Bin - + — ^,1 
1 2 2 ' 



> 



d + l 



>l-d^ 1 ), 
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and therefore < 1/d. This proves the claim. ■ 

Proof of Lemma 4.5: As we pointed out earlier, although the subtrees (T t **)*il^ +1 are correlated by the construc- 
tion of the islands, they are independent of the substitution process. By forcing (randomly) the subtrees (T t **)*,"t!^ + x 
to be (d — 2)-ary and fixing the adversarial nodes to 1, we restore the i.i.d. nature of the sites, from which the result 
follows. ■ 

Proof of Lemma 4.7: This follows from Lemma 4.5, the independence of A and 0, and three applications of 
Hoeff ding's lemma. ■ 

Proof of Lemma 4.8: Note that 

1 m ^ m 

Corr(A,y) = -J>i>(fc) = - J><> (l/i> <*> < 

rn ^ — * rn ^ — * 

Hence, 



m * — ' ' m 

i=l i=l 



in in 

| Corr(A,y) - Corr(A, Y)\ < - V(l - = 1 - - Y>i} ( Wi ). 

t=i i=i 

Now notice by case analysis that 

(Zi)( Wi ) > (AiX^i) - l(A i >=-l " %>=-!• 

This proves the first claim. The second claim follows from the bounds in Lemma 4.7. ■ 
Proof of Lemma 5.1: For the first claim, note that 

E[Corr(<,0] = 9 2 s = (1 - 2 Ps ) 2 , 

where we used that 1) there is no indel in the sites [t + 1 : t + a] between v and 2) that the sites are perfectly 
aligned; and 3) that the substitution process is independent of the indel process. We also used the well-known fact 
that the # s 's are multiplicative along a path under our model of substitution [SS03]. The result then follows from 
Hoeffding's inequality. 

For the second claim, because the anchors are now misaligned the t'-th term in Corr( I / r ?, ^) for t' € [t + 1 : 
t + a] is the variable {x\, +s .( r ))(a^, +s ( r )_i) which is uniform in {—1, +1}- In particular, we now have 

E[Corr(y;,^')] = 0. 

The result follows from the method of bounded differences with respect to the independent vectors 



yt+a 
=f 



Proof of Lemma 5.2: This follows from Lemmas 4.8 and 5.1 and the triangle inequality. ■ 

D Completing the Proof of the Main Theorem 

Having conditioned on the event £, we justify the correctness of our reconstruction method via the following 
induction. The top level of the induction establishes Theorem 1.3. 

Induction hypothesis. Consider a parent v in T* ; in particular, v is stable. We assume that the following conditions, 
denoted by (*), are satisfied: For all children i G [d] of v belonging to T* 



20 



1. Alignment. For all children %' of i with i' G T* and all r = 1, . . . , kjt — 1, 



Mr) = Si'(r). 



(2) 



(This condition is trivially satisfied for values of rl that are larger than the sequence length of i'.) 



2. Reconstruction. Moreover, we have ki = ki and for all t = 1, . . . , ki, the following holds: 



Let Li be the leaves below i with = |Lj|. Let H be the level of v. Let L\* be the gateway 
leaves for site t. For u G L^* let F u (t) be the position of site t in u. Note that x\ can be written 
as x\ = Maj^ _1 (zi, . . . , z ni ), where zj is either ft or for an appropriate function \>j. Our 
hypothesis is that 




In particular, the ancestral reconstruction Xi dominates the adversarial reconstruction X[. 
The base case where v is a leaf is trivially satisfied. 

Alignment. We begin with the correctness of the alignment. 

Lemma D.l (Induction: Alignment) Assuming £ and (*), the algorithm infers Si correctly for all children i G [d] 
which are also in T*, that is, (2) holds for v. 

Proof of Lemma D.l: Let II denote the set of children of v in T*. The proof follows by induction on r. The base 
case r = is trivial. Assume correctness for r — 1. 

If there is no indel in any of the children i G II between the sites (r — l)£ and rl of v, then under £, (*) and 
Lemma 5.2 we have II C G r . In that case, for all i G II we have Si(r) = §i(r — 1) = Si(r — 1) = Sj(r), where the 
second equality is from (*). 

If there is an indel operation in island r, then since v is stable only one indel operation occurred in one child. 
Denote the child with an indel by j. Assume the indel is a deletion. (The case of the insertion is handled similarly.) 
If j is not in T* we are back to the previous case. So assume j is in T*. Again, from £, (*) and Lemma 5.2 the 
other children in T* are added to the set G r , and the shift value will be computed correctly for them. Moreover by 
(*), for every i G II — {j}, 



fi(r£ + l) 



rl + 
rl + 
rl + 



1 + Si (r) 
1 + Si(r) 
1 + Si(r - 



1) 



which is the starting point of A l r . Also, 



fj(rl+l) 



rl+1 + Sj(r) 
rl+ 1 + sj(r - 
rl + 1 + §j (r - 
rl + Sj(r- 1), 



1)-1 
1)-1 



which is the stalling point of D 3 r . Thus according to Lemma 5.2 D 3 r matches A\ for alH G IT fl G r . As there are 
d — 2 children in II n G r , we get that the algorithm sets 



Sj(r) = sj(r — 1) — 1 = Sj(r — 1) — 1 = Sj(r) 



as required. Note also that in this case, according to Lemma 5.2 again, A 3 r does not have high correlation with A l T 
for any i G II fl G r , and thus we will consider P r and D J r . Similarly, P r does not have high correlation with A\, for 
any i G II n G r , and thus we will not try to set Sj(r) twice. ■ 
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Ancestral reconstruction. We use Lemma D. 1 to prove that the ancestral reconstruction dominates the adversarial 
reconstruction. In the algorithm, we perform a sitewise majority vote over the children of v in G r (these are the 
aligned children — see the description of the algorithm in Figure 1). For notational convenience, we assume that in 
fact we perform a majority vote over all children but we replace the states of the children outside G r with (J. 

Lemma D.2 (Induction: Reconstruction) Assuming £, (*) and the conclusion of Lemma D.l, (3) holds for v. In 
particular, the ancestral reconstruction X v dominates the adversarial reconstruction X' v . 

Proof of Lemma D.2: The second claim follows from the first one together with the construction of the adversarial 
process and the mono tonicity of majority. 

As for the first claim, by Lemma D.l for each site of v there are d — 2 uncorrupted children islands containing 
this site such that the children are also in T*. In particular, the d—2 corresponding sites in the children are correctly 
aligned. Moreover, by the induction hypothesis, each corresponding site in the children satisfy (3). By taking a 
majority vote over these sites we get (3) for v as well. 

A small technical detail is handling the case where the last island has less than a sites, and thus does not contain 
an anchor. However, in this case, if the father is stable then there are no indel operations at all in the last island, and 
therefore aligning it according to the previous one gives the right result. ■ 
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